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Abstract 

Potassium-competitive acid blocl<ers (P-CABs) are liiglily safe and active drugs targeting H"^,K'^-ATPase to cure acid-related 
gastric diseases. In this study, we for the first time investigate the interaction mechanism between the protonated form of P- 
CABs and human H^,K^-ATPase using homology modeling, molecular docking, molecular dynamics and binding free energy 
calculation methods. The results explain why P-CABs have higher activities with higher pKa values or at lower pH. With 
positive charge, the protonated forms of P-CABs have more competitive advantage to block potassium ion into luminal 
channel and to bind with H^,K^-ATPase via electrostatic interactions. The binding affinity of the protonated form is more 
favorable than that of the neutral P-CABs. In particular, Aspl39 should be a very important binding site for the protonated 
form of P-CABs through hydrogen bonds and electrostatic interactions. These findings could promote the rational design of 
novel P-CABs. 
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introduction 

The gastric H"'",K^-ATPase (proton pump), primarily responsi- 
ble for gastric acid secretion, is the key therapeutic target for the 
ulcer diseases such as duodenal ulcers, gastric ulcers, gastro 
esophageal reflux disease (GERD), ZoUinger-EUison syndrome 
(Z-E), and gastritis [1-3]. As a member of the P2-type ATPase 
family, H^,K^-ATPase is a dimeric heterodimer composed of a 
subunit of about 1035 amino acids with 10 transmembrane (TM) 
segments and /J-subunit glycoprotein with 290 amino acids [4,5]. 
By cyclic phosphorylation and dephosphorylation of the catalytic 
subunit, H''",K''"-ATPase undergoes conformational changes 
between Ei and E2. With phosphorylation, H^,K^-ATPase Ei 
conformation binding hydronium ion EiP'H^O'*' changes to 
E2P'H30^ form. After release of H^O"*" and binding of K""" on 
the extracytoplasmic surface, the E2P»K^ conformation is formed 
and then converts to the EjK conformation with the dephosphor- 
ylation. The EiK conformation releases K""" to the cytoplasmic 
side, then rebinding of H^O^ occurs to complete the transition 
cycle [6]. The H+, K+- ATPase engages in 2K""/2H^/1ATP 
electroneutral ion exchange, generating a million-fold H''"-gradient 
across the mammalian canalicular membrane of the parietal cell 
[7,8]. 

H''",K''"-ATPase inhibitors include two classes: proton pump 
inhibitors (PPIs) and potassium-competitive acid blockers (P- 
C ABs) [9,10]. PPIs such as omeprazole, lansoprazole, rabeprazole, 
pantoprazole, tenatoprazole and leminoprazole, are irreversible 
inhibitors of H''",K''"- ATPase, which form a covalent complex with 
the protein at specific cysteine residues [11,12]. Although currently 
recognized as the most effective drugs for the treatment of acid- 
related diseases, PPIs exhibit a delayed onset of acute effect and 



achieve full effect only slowly and incrementally over several dose 
cycles [13]. Therefore, many patients with GERD symptom are 
unsatisfied with PPIs treatment [14,15]. An alternative to PPIs is 
P-CABs, which reversibly inhibit gastric H''',K'''-ATPase by 
competing with the K"*" on the luminal surface [1]. After oral 
dosing, P-CABs rapidly achieve high plasma concentrations with a 
fast onset of action [13]. 

Now several P-CABs including SCH28080 [16-18], Sorapra- 
zan [19,20], Revaprazan [21-23], AZD0865 [24,25] and TAK- 
438 [26-29] (Figure 1) were developed by pharmaceutical 
companies. Among them Revaprazan was used clinically in 
2007 for the treatment of duodenal ulcer, gastric ulcer and 
gastritis, and is undergoing phase III clinical studies for the 
treatment of GERD. AU P-CABs are weak bases. SCH28080, 
AZD0865 and TAK-438 have pKa values of 5.6 [30], 6.1 [24] 
and 9.37 [27], respectively (Table 1). Gedda et al. [24] reported 
that at pH 7.4 AZD0865 concentration-dependendy inhibited 
K"^-stimulated H'^,K"^-ATPase activity (IC50 = 1.0 nM) but was 
more potent at pH 6.4 (IC5q = 0.13 |xM). The theoretical proton- 
ated AZD0865 is approximately 33% at pH 6.4 and less than 5% 
at pH 7.4. SCH28080 has also been reported to be weaker under 
neutral conditions (IC5o = 0.14 |J,M at pH 6.5; IC5o = 2.5 \lM at 
pH 7.5). In comparison, because of its high pKa value, TAK-438 
should be protonated instantly and exert a potent inhibitory 
activity even in a neutral environment (IC50 =^ 0.019 |lM at 
pH 6.5; IC5„ = 0.028 at pH 7.5) [27]. According to the pKa 
calculation using ACD/I-Lab [31], SCH28080 and TAK-438 are 
8.08% and 98.36% protonated at pH 6.5, while 0.87% and 
94.67% protonated at pH 7.5 (Table 1). Thus, the protonated 
form would be the highly active form of P-CABs. But the 
mechanism of interaction between protonated P-CABs and 
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(AstraZeneca) 



AZD0865-1 



Figure 1. Chemical structures of P-CABs with their protonated form. 

doi:1 0.1 371 /journal.pone.0097688.g001 



H^jK^-ATPase is not known in detail, which hinders the 
development of novel P-CABs. 

So far the structural investigations of H^,K"'"-ATPase have 
lagged behind the pharmacological studies. The structure of 



human gastric H^,K^-ATPase is unknown, and the structure of pig 
gastric H"'^,K^-ATPase is poorly defined, being currently limited to 
a resolution of 7 A (PDB code: 3IXZ [32], resolution: 6.5 A; PDB 
code: 2XZB [33], resolution: 7 A). So the aim of this study is to 



Table 1. IC50, pKa (reference) and pKa (calculation) values of P-CABs. 





P-CABs 






pKa(Ref.) 


pKa(Cal.) 


protonated form (+1) percentage/%° 




pH<7.0 


pH£7.0 






pH 6.5 


pH 7.5 


SCH28080 


0.14 (pH 6.5) 


2.5 (pH 7.5) 


5.6[ref.30] 


5.85±0.10 


8.08 


0.87 


Soraprazan 




0.1 (pH 7.0) 




7.11 ±0.70 


95.50'' 


72.84 


Revaprazan 


0.35 (pH 6.1) 






7.26±0.10 


69.89' 


10.42 


AZD0865 


0.13 (pH 6.4) 


1 .0 (pH 7.4) 


6.1[ref.24] 


6.50±0.10 


33.00 


5.00 


TAK-438 


0.019 (pH 6.5) 


0.028 (pH 7.5) 


9.37[ref.27] 


9.06±0.10 


98.36 


94.67 



^protonated form (-t-1) percentages were calculated using ACD/I-Lab except for AZD0865 which from reference [24]; b. protonated form percentage of Soraprazan is 
87.07% at pH 7.0; c. protonated form percentage of revaprazan is 69.89% at pH 6.1. 
doi:l 0.1 371/journal.pone.0097688.t001 
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model H^jK^-ATPase structure by homology modeling and for 
the first time to investigate the interactions between the protonated 
form of P-CABs and H"*^,K*-ATPase using molecular docking, 
molecular dynamics and MM/GBSA calculation methods. 

Materials and Methods 

Homology modeling 

The sequences of the human (1035 amino acids) and pig (1033 
amino acids) gastric H"'",K"'"-ATPase receptor were taken from the 
Swiss-Prot Database (ID: P20648 and P09626) [34]. Through 
BLAST online method (http:/ / blast.ncbi.nkn.nih.gov/Blast.cgi) 
[35] from the Protein Data Bank [36], the crystal structure of 
Na'^,K"*'-ATPase in the E2P state (PDB code: 2ZXE) [37] was used 
as a template. The sequence alignment was performed with the 
ClustalW2 algorithm [38]. The homology model of human 
H"',K"'-ATPase was generated using MODELLER9v4 [39]. The 
resultant structure of the H"'",K"'"-ATPase was subject to the Protein 
Preparation Wizard module in Schrodinger [40] as follows: adding 
hydrogens, assigning partial charges, and minimizing using the 
OPLS-2005 force field [41] until RMSD 0.30 A. The final 
optimized model was validated using the program PROCHECK 
[42] to assess the quahty of the stereochemistry of the protein 
structure. 




Figure 3. The human gastric H'^,K'^-ATPase model. 

doi:10.1371/journal.pone.0097688.g003 
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Figure 2. Sequence alignment results between human H^,K^-ATPase and the template Na%K^-ATPase (2ZXE). The residues with 
identical, strong and weal< similarities in K*-ATPase are shown in dark blue, blue and light blue color background, respectively. The alpha helical, 
sheet and coil of secondary structure in the template are colored by pink, light purple and light green. 
doi:1 0.1 371/journal.pone.0097688.g002 
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Preparation of ligands 

LigPrep of Schrodinger software suit [43] was used for the 
preparation of ligands: generating 3D structures from 2D (SDF) 
representation, and performing a geometry minimization. Tlie 
ligands were subjected to energy minimization using MacroModel 
module of Schrodinger with Merck Molecular Force Field 
(MMFFs). Truncated Newton Conjugate Gradient (TNCG) 
minimization method was used with 500 iterations and conver- 
gence threshold of 0.05 (kj/mol). While Epik [44] was used to 
generate possible ionization states at pH 7.0±1.0. 

Molecular docking 

The docking simulations were performed using the software 
Glide pCP mode) [45,46]. Previous biochemical and mutagenesis 
studies [18,33,47-50] suggest that Cys813 in pig H'^,K'^-ATPase 
(corresponding to Cys815 in human H''",K"''-ATPase) is the key 
amino acid residue in the luminal cavity. Therefore, dimensions 
for the cubic boundary box centered on the centroid of Cys8 1 5 
were set to 20 A x 20 A x 20 A. The scaling factor for protein van 
der Waals radii was 1.0 in the receptor grid generation. After 
docking calculations, for each ligand, the best pose was chosen and 
scored using the proprietary GhdeScore function. 

QM/MM optimization 

The final docking complexes were energetically optimized by 
QM/MM method. QM/MM calculations were carried out using 
the QSite program [51,52] of the Schrodinger suite. The ligands 
were defined as QM region calculated by the density functional 
theory DFT/B3LYP (6-3 IG* basis set). The receptor as MM 
region was minimized with Truncated Newton algorithm (max- 
imum cycles as 1000; gradient criterion as 0.01). The OPLS 2005 
all-atom force field was employed. 

Molecular dynamics 

The optimized docking models were subjected to molecular 
dynamics simulations using Desmond [53,54]. The system was 
embedded in a fuUy hydrated POPC (l-palmitoyl-2oleoyl-OT- 
glycero-3-phosphatidylchlorine) bilayer membrane and solvated 
with an orthorhombic box of SPG water molecules (buffer 
distance: 10 Ax 10 Ax 12 A). Counter-ions (Na"^ were added to 



neutralize the system and a physiological salt concentration of 
0.15 M NaCl was introduced. The final system was composed of 
approximately 130,000 atoms. Each model was equilibrated in 
MD for 10 ns. Then two K"*" were added in the luminal domain 
near channel (close to Asp 134) replacing two Na"*", and potassium 
ion competition molecular dynamics simulations were carried out 
for 10 ns. Before the simulation, the models were relaxed as 
follows: (1) two minimization steps (restraining the solute and 
unrestrained minimization) with maximum runs of 2000 and the 
convergence threshold for minimization set to 1 kcal/ mol/ A. The 
minimization method was a hybrid of the steepest decent and 
limited-memory Broyden-Fletcher-Goldfarb-Shanno (LBFGS) al- 
gorithms; (2) after minimization, the simulation in the NVT 
ensemble was run restraining all solute heavy atoms with 
temperature of 10 K for 20 ps, using Berendsen thermostat; (3) 
a simulation in the NPT ensemble restraining all solute heavy 
atoms with temperature of 10 K and 300K for 50 ps, respectively; 
(4) a simulation in the NPT ensemble, no restraints, with 
temperature of 300 K and simulation time of 200 ps. Then 
10 ns MD production runs (time step: 2.0 fs) were performed 
through NPT ensemble at 300 K with 1.0132 bar pressure. 
Smooth particle mesh Ewald method (Ewald tolerance: le-09) was 
employed to treat the long-range electrostatic interactions and a 
9 A radius cut off was used for coulombic short range interactions. 
The energies and frames of each trajectory were recorded every 1 
ps and 5 ps, respectively. After 1 0 ns MD simulations with two K"^, 
SCH28080 and TAK438-1 systems with H'",K'"-ATPase and 
POPC membrane in 3M KCl resolution were equihbrated in MD 
for 40 ns and then run 100 ns disassociation molecular dynamics, 
respectively. MD trajectory analysis was performed using Des- 
mond utilities and VMD [55]. The ligand-protein complexes were 
visualized using PyMOL [56] and analyzed with Ligand 
Interactions module embedded in Maestro 9.3 [57]. 

MM/GBSA calculations 

For each system, binding free energy (AGh;„a) calculations were 
performed for 100 snapshots extracted from the last 1 ns stable 
MD trajectory using molecular mechanics-generalized Born 
surface area (MM/GBSA) method. MM/GBSA procedure in 
Prime program [58,59] was used to calculate AGi,ind of the docked 
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Figure 4. RMSDs for the backbone atoms of the complexes: (A) P-CABs; (B) protonated form of P-CABs. 

doi:1 0.1 371 /journal.pone.0097688.g004 
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Figure 5. Time evolutions of the distances between potassium ion (K1) and nitrogen atom of P-CABs (N 15475 of SCH28080 or 
N15494of TAK-438-1) (A). Tine distances beyond 35 A of TAK-438-1 are not shown in plot. The interaction modes of SCH28080 system at 0 ns (B) 
and 7.570 ns (C), and TAK-438-1 system at 0 ns (D) and 3.180 ns (E), with the potassium ion represented by blue ball. 
doi:1 0.1 371 /journal.pone.0097688.g005 



inhibitors according to the following equations [60] : 

^ ^biiid = ^EmM + a Gjo/v ( 1 ) 



AGw = A£-MM + AG,w.-rAS (2) 

Where AEmm is the difTerence of the gas phase MM energy 
between the complex and the sum of the energies of the protein 
and inhibitor, and includes AEii-„en,aj (bond, angle, and dihedral 
energies), AEEied (electrostatic), and AEvdw (van der Waals) 
energies. AG,„iv is the change of the solvation free energy upon 



binding, and includes the electrostatic solvation free energy AGgb 
(polar contribution calculated using generalized Born model), and 
the nonelectrostatic solvation component AGg^ (nonpolar contri- 
bution estimated by solvent accessible surface area). TAS is the 
change of the conformational entropy upon binding. This term 
was calculated using normal-mode analysis Rigid Rotor Harmonic 
Oscillator (RRHO) contained in MacroModel module [61]. 
AG neglects the effect of entropy contributions, while AGynd 
includes contributions from loss of ligand translational, rotational 
and vibrational entropy (TAS). 
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Figure 6. RMSF of each residue in SCH28080 and TAK-438-1 complexes. 
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Results and Discussion 

H^K^-ATPase homology model 

The three-dimen.sional structure of Na^,K^-ATPase in the E2P 
state (PDB code: 2ZXE; resolution: 2.4 A) [37] was selected as a 
template, which share 63.9% identity and 78.7% similarity to 
human H''",K''"-ATPase on the basis of sequence alignment analysis 
(Figure 2). The human gastric H''',K"''-ATPase model is shown in 
Figure 3. It comprises ten transmembrane helices (TMI-TMIO), 
in which binding sites are located, and three cytoplasmic domains: 
the nucleotide-binding (N), phosphorylation (P) and actuator (A) 
domains [33]. The stereochemistry of the homology model was 
assessed using Ramachandran plot, which indicates that 93.8% of 
the residues were located in the most favored zones, 6.1% in 
allowed regions, 0.0% in generously allowed regions and 0.1% in 
disallowed regions (Figure SI). The dihedrals, covalent and overall 
G-factors of this model are -0.17, 0.01 and -0.08, respectively. The 
PROCHECK G-factor should be above -0.5 ideally for the 
homology model to have a good acceptance score. 



Molecular dynamics simulations 

After molecular docking and QM/MM optimization (Glide 
docking Gscore, Glide energy and QM/MM energy were listed in 
Table SI), MD simulations for ten P-CAB-H+,K"^-ATPase 
complexes with two K""" in NaCl aqueous solution were run for a 
duration of 10 ns. To explore the dynamic stability of complexes 
and to ensure the rationality of the sampling method, root-mean- 
square deviations (RMSD) for the backbone atoms from the 
starting structure were analyzed, as shown in Figure 4. As can be 
seen in the plots, after 8 ns, the RMSD of each system tends to 
converge, indicating that the systems are stable and equihbrated. 
During the whole MD simulation process, the potassium ions were 
blocked into the luminal channel by all P-CABs, which acted just 
like goalkeepers. The distances between one of potassium ions (Kl) 
and the nitrogen atom (N 15475) of -ON in SCH28080 (or 
nitrogen atom (N15494) of-NH2"^CH3 in TAK-438-1) are shown 
in Figure 5 with the interaction modes at 0 ns and the time of the 
shortest distances (SCH28080: 7.570 ns, distance 2.744 A; TAK- 
438-1: 3.180 n.s, distance 6.218 A). Due to the electrostatic 
repulsive interaction, the protonated form (H-1 charged) can block 
K""" more easily than the neutral P-CABs. From 4.218 ns, the 



Table 2. RMSF values (A) of important amino acid residues in different P-CABs system. 







Leu 143 


Ala337 


TyrSOl 


Leu81 1 


Cys815 


SCH28080 


1.267 


1.198 


1.038 


1.544 


1.196 


Soraprazan 


1.205 


1.138 


0.998 


1.280 


1.243 


Revaprazan 


2.403 


1.489 


1.483 


1.905 


1.212 


AZD0865 


1.075 


1.173 


1.069 


1.286 


1.239 


TAK-438 


1.109 


0.985 


0.837 


1.081 


0.981 


SCH28080-1 


1.831 


1.913 


2.174 


2.134 


1.758 


Soraprazan-1 


1.637 


1.260 


1.005 


1.193 


1.340 


Revaprazan-1 


1.391 


1.025 


0.867 


1.226 


0.964 


AZD0865-1 


1.338 


1.197 


0.979 


1.544 


1.062 


TAK-438-1 


1.642 


0.991 


0.949 


1.067 


0.865 



The smallest RMSF values of residues among P-CABs complexes are bold. 
doi:l 0.1 371/journal.pone.C097688.t002 
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Charged (negative) j Metal — n-catjon 

Charged (positive) hzO Water ~ H-bond (backbone) 
Polar J Hydration site - H-bond (side chain) 

Hydrophobic ^ Displaced hydration site ^ Metal coordination 

Glycine ^ n-n stacking Solvent exposure 



Figure 7. Interaction modes of P-CABs witKi IH^K^-ATPase. (Aj 

SCH28080; (B) SCH28080-1; (C) TAK-438; (D) TAK-438-1; (E) Soraprazan; 
(F) Soraprazan-1; (G) Revaprazan; (H) Revaprazan-1; (I) AZD0865; (J) 
AZD0865-1 

doi:10.1371/journal.pone.0097688.g007 

distance between K*** and TAK-438-1 is larger than 20 A, while 
this time point is 9.421 ns for SCH28080 system (Figure 5A). 
Furthermore, analyses of root-mean-square fluctuation (RMSF) 
versus the residue number for SCH28080 and TAK-438- 
1 complexes are illustrated in Figure 6. All complexes share similar 
RMSF distributions and similar trends of dynamic features. They 
are relatively rigid in the active site region (residues Leu 143 in 
TM2, Ala337 in TM4, Tyr801 in TM5, Leu811 in the TM5-6 
loop, and Cys815 in TM6) as reported in the literatures 
(residue IDs in pig H ,K -ATPase correspondingly are Leu 141, 
Ala335, Tyr799, Leu809, and Cys813). The fluctuation of the 
active site for TAK-438-1 complex (except for Leu 143) is smaller 
than SCH28080 complex. RMSF values of the important amino 
acids in different P-CAB system are listed in Table 2. 

Interaction modes of P-CABs with H^K'^-ATPase 

To investigate P-CAB interactions in the binding site, the 
average structures from last 1 ns MD trajectory were compared 
(Figure 7). After MD simulation, the pose of SCH28080 was very 
similar to that reported by Abe et al. [33] . The binding site is in the 
luminal cavity, surrounded by Ala337, Tyr801 and Cys815, which 
have hydrophobic interactions with SCH28080 (Figure 7A). The 
hydrogen bond (H-bond) interactions, which play an important 
role in P-CABs binding to H'^,K"^-ATPase, are hsted in Table 3 
and shown in Figure 7. Oxygen atom of SCH28080 and nitrogen 
atom of AZD0865 formed hydrogen bonds with Tyr930. TAK- 
438 was found to form H-bond with Tyr801, and has Tt- 7i stacking 
interaction with Phe334 and 7i-cation interaction with Arg330. 
The protonated forms have different interaction modes compared 
to the neutral P-CABs. There are strong H-bond interactions 
between Aspl39 and TAK-438-1 (distance 1.727±0.012 A) and 
Soraprazan-1 (distance 1.623±0.006 A and 2.133±0.015 A). 
Asp 139 also has negative charged interactions with aU protonated 
P-CABs (-Hi charged). AZD0865-1 interacts with Aspl38 via two 
H-bonds (distance 1.695±0.009 A and 1.903±0.013 A). Both 
Soraprazan-1 and Revaprazan-1 have H-bonds with Phe990. 
TAK-438-1 formed hydrogen bonds with Asp 139, Tyr927 and 
Asn991 simultaneously. 

Binding free energy of P-CABs 

The binding free energies of aU ten systems were calculated by 
MM/GBSA method. As hsted in Table 4, the AG bmd and AGbind 
values all show that the binding afiinity of protonated P-CABs is 
more favorable than that of neutral P-CABs. The order of favorable 
binding interaction is TAK-438- l>Soraprazan-l>Revaprazan-l> 
SCH28080-1>AZD0865-1>TAK-438>AZD0865>SCH28080« 
Soraprazan>Revaprazan, which is consistent with the experimental 
results at different pH. According to the protonated form percentage 
of P-CABs in Table 1, the binding free energies (AG|,„jd) of 
protonated and neutral mixtures of P-CABs at different pH values 
were calculated. The correlation coefficient R between AGi,;,,^ of the 
mixture and pICso (negative logarithms of 50% inhibition concen- 
tration) is -0.837 (Figure 8). 

To estimate which energy term has most impact on the binding 
affinities, the four individual energy components (AEEkd, '^^EvdWj 
AGgb, and AG,sa) were carefully compared. From Table 4, it can 
be seen that both the van der Waals (AEvdw) and the electrostatic 
(AEEiect) contributions are essential for P-CABs binding to H , 
K''"-ATPase. For the neutral form, the contributions of AEvdw are 
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Table 3. Hydrogen bond analysis after MD simulation. 



P-CABs 


Atoms 


Residues 


Atoms 


Distance (A)* 


SCH28080 


01 


Tyr930 


HH 


1.986±0.016 


TAK-438 


02 


Tyr801 


HH 


1.947±0.013 


TAK-438-1 


H16 


Aspl39 


0D2 


1.727±0.012 




H17 


Tyr927 


OH 


2.188±0.019 




N2 


Asn991 


HD22 


2.284±0.020 


Soraprazan 


H20 


Aspl39 


0D2 


1.767±0.009 


Soraprazan-1 


HI 2 


Asp! 39 


0D1 


2.. \ JJ_U.U 1 D 




H26 


Aspl39 


ODl 


1.623 ±0.006 




H20 


Phe990 


0 


1.976±0.014 


Revaprazan-1 


H24 


Phe990 


0 


1.898±0.012 


AZD0865 


H14 


Glu902 


0E2 


1.931 ±0.025 




N2 


Tyr930 


HH 


2.118±0.018 


AZD0865-1 


HIS 


Aspl38 


ODl 


1.903±0.013 




H27 


Aspl38 


ODl 


1.695±0.009 


*The average distance with standard error (SE = 


standard deviation/N^^^} of hydrogen bond 


in the last 1 ns MD. 





doi:l 0.1 371/journal.pone.0097688.t003 



more favorable than AEgiect term. But for protonated P-CABs, the 
major favorable contributor Ls AEeicci term, which is far larger 
than for the neutral form. Although polar solvation (AGgb) term of 
protonated form is also far larger than for the neutral form, which 
could offset AEEiect, the net electrostatic contributions (AEEif,tt+ 
AGgb) are more favorable for the protonated compounds. Thus, 
the electrostatic contribution is more crucial for distinguishing the 
binding afTmities among P-CAB complexes. Except for Sorapra- 
zan-1 and AZD0865-1, nonpolar solvation term (AGg^) contribute 
favorably to binding. Among all P-CABs, AEeicci and AGgA of 
TAK-438-1 is the most favorable (- 145. 88± 1.80 and 
— 9.61±0.30 kcal/mol). In addition, the contributions of the 
conformational entropy (TAS) are very similar for the protonated 
compounds (-2.34±0.21 3.32±0.29 kcal/mol). 



8.5-, 




(PH 7.5) 

1 — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' 
-70 -65 -60 -55 -50 -45 -40 -35 -30 

Binding free energy (kcal/mol) 

Figure 8. The relationship between binding free energy and 
pICso of P-CABs at different pH. 

doi:1 0.1 371/journal.pone.0097688.g008 



The binding free energy between the protein and P-CABs was 
decomposed into the contribution of each residue, which provides 
quantitative information of the key residues related to the detailed 
binding mechanism. From the energy comparison of residues in 
binding sites (Figure 9), it can be observed that there is the distinct 
difference between protonated and neutral forms of P-CABs. The 
energies of residues Tyrl42, Tyr804, Leu811, Leu813, Gly814, 
and Tyr930 interacting with SCH28080-1 (-0.40±0.10, -0.39± 
0.15, -0.31±0.02, -0.55±0.07, -0.33±0.12, and -0.20± 
0.04 kcal/mol) are less favorable than those of SCH28080 
complex (-1.79±0.07, -3.18±0.34, -3.46±0.13, -2.26±0.10, 
-5.35±0.09, and -3.02±0.49 kcal/mol). But the interactions 
between SCH28080-1 and residues Aspl39, Asnl40, Leul43, 
Val333, Met336, Ala337, Val340, Tyr801, Cys815, and Ile816 
of H^K-'-ATPase (-2.92±0.62, -3.96±0.45, -5.46±0.25, 
-2.35±0.20, -5.55±0.48, -3.31±0.42, -3.09±0.17, -5.94± 
0.35, — 6.56±0.31, and — 5.70±0.15 kcal/mol) are more favor- 
able than those of SCH28080 (-0.16±0.24, 0.08±0.03, -2.73± 
0.08, -0.95±0.08, -0.60±0.05, -0.95±0.08, -0.07±0.02, -1.89± 
0.13, -3.07±0.22, and -2.49±0.08 kcal/mol) (Figure 9A), which in 
particular include the key residues Leu 143, Tyr801, and Cys815 (all 
above -5.00 kcal/mol for SCH28080-1). Tlris might be the reason 
tliat the activity of SCH28080 at pH 6.5 is dramatically higher than 
that at pH 7.5. Due to strong hydrogen bond and electrostatic 
interactions, the energy contributions of Asp 1 39 to protonated form are 
all more favorable than those to neutral P-CABs (— 14.63±0.51 kcal/ 
mol to Soraprazan-1, -1 1.29±0.35 kcal/mol to Revaprazan-1, 
— 9.57±0.26 kcal/mol to AZD0865-1). In particular, the highest value 
(-23.28±0.64 kcal/mol) is reached for TAK-438-1. Even tliough less 
favorable than their protonated forms, the interactions with Asp 1 39 are 
also strong to Soraprazan (— 5.69±0.77 kcal/mol) by H-bond and to 
AZD0865 (-7.14±0.41 kcal/mol) by electrostatic interactions. 
Asp 139 could be the key residue in the binding sites of P-CABs, 
which was not paid attention to previously. TAK-438 has more 
favorable interactions with the key residues Val333, Tyr801, 
Leu811, and Cys815 (-4.88±0.13, -3.59±0.16, -2.86±0.09, 
— 2.84±0.17 kcal/mol). However, tlie binding free energy of TAK- 
438-1 is higher than TAK-438 because TAK-438-1 has strong 
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interactions witli Asp 139 (H-bond and electrostatic interactions), 
Asnl40 (polar interaction, -8.13±0.25 kcal/mol), Tyr927 (H-bond 
and hydrophobic interactions, — 11.51±0.56 kcal/mol), and Asn991 
(H-bond and polar interactions, — 12.85±0.53 kcal/mol) (Figure 7D). 
Soraprazan-1 and Revaprazan-1 have almost the same features 
compared to TAK-438-1. Compared to their neutral form, they have 
more favorable interactions with residues Aspl39, Asnl40, Tyrl42, 
Leul43, Tyr927, Tyr930, Phe990, Asn991, and Phe992. hi addition, 
there is strong hydrophobic interaction between Soraprazan-1 and 
LeuBll (-5.72±0.24 kcal/mol). Different to tlie other P-CABs, 
Revaprazan-1 has hydrophobic interactions with Cysl22 (— 1.35± 
0.08 kcal/mol), Aal25 (-2.15±0.21 kcal/mol), and electrostatic 
interaction with Glnl29 (-3.02±0.34 kcal/mol). As well as hydro- 
phobic interactions with Pro810 (-2.97±0. 13 kcal/mol), LeuSll 
(-3.26±0.24 kcal/mol), and Cys815 (-4.22±0.26 kcal/mol), 
AZD0865 interacts with Glu902 (-9.04±0.70 kcal/mol) and 
Tyr930 (-8.89±0.37 kcal/mol) by H-bonds, with Gln926 
(— 5.53±0.23 kcal/mol) by polar interaction, and with Apl39 by 
electrostatic interaction. Positively charged Arg330 has unfavorable 
interactions with AZD0865 and Soraprazan-1 (4.37±0.45 and 
2.19±0.15 kcal/mol). Because of the strong H-bond and electrostatic 
interactions with nitrogen atoms, the energy contribution of Apl38 is 
high (-15.82±0.41 kcal/mol) for AZD0865-1, second only to the 
contribution of Ap 1 39 for TAK-438- 1 . The results demonstrate that 
the hydrogen bond and electrostatic interactions with Apl39 
(or Asp 138) should be very important for P-CABs binding to H^, 
K*-ATPase, which are in agreement with tlie reference [62] . Using the 
competitive inhibitor 8-[(4-azidophenyl) methoxy]-l-trithiomethyl-2,3- 
dimethylimidazo-(l,2-a) pyrimidium iodide (mDAZIP, a photoactiva- 
ble compound derived from SCH 28080), Munson et al. [62] suggested 
the binding site was on the luminal side between Glnl29 and An 140 
(between Gkil27 and Anl38 in the TMl-2 loop of pig H"*", 
K-^-ATPase). 
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Disassociation molecular dynamics 

After 100 ns molecular dynamics (RMSD for the backbone 
atoms and RMSF for the residues of SCH28080 and TAK-438- 
1 complexes are illustrated in Figure S2 and S3), SCH28080 
occurred partial disassociation from H''",K''"-ATPase and AG|,i,j;j of 
SCH28080 decreased to -27.88±0.74 kcal/mol (Table 4). 
Because water molecules have went into the ion channel, the 
interactions between SCH28080 and residues Val333 to Val340 in 
TM4 were cut off (Figure lOA and Figure 11). In contrast to 
SCH28080, AGi,i,jd of TAK-438-1 after 100 ns MD was - 66.48 ± 
1.87 kcal/mol (Table 4) and similar to that in 10 ns MD. TAK- 
438-1 stm interacted with Aspl39 (-25.48±0.90 kcal/mol), 
Tyr927 (- 1 1.81 ± 1.28 kcal/mol) and Asn991 (-9.33± 
0.31 kcal/mol) by H-bonds and blocked water molecules into 
the ion channel (Figure lOB and Figure 1 1). Therefore TAK-438- 
1 could have the long dwell time and disassociate from H''", 
K*-ATPase very slowly, which make it a likely competitor for 
PPIs. 

Conclusions 
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In the present study, to clarify the interaction mechanism 
between protonated P-CABs and H"'",K"''-ATPase, MD simulations 
and MM/ GBSA binding free energy calculations were conducted. 
There is a distinct difference in interaction mode between 
protonated and neutral P-CABs. With positive charges on nitrogen 
atoms, protonated forms of P-CABs have a competitive advantage 
in blocking potassium ion into the luminal channel and binding to 
H^,K*-ATPase mostly via electrostatic interactions. According to 
the binding free energy analysis, the binding affinity of protonated 
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Figure 9. The comparison of energy decomposition for residues in binding sites of P-CABs. (A) SCI-128080 and SCI-128080-1; (B) TAK-438 
and TAK-438-1; (C) Soraprazan and Soraprazan-1; (D) Revaprazan and Revaprazan-1; (E) AZD0865 and AZD0865-1 
doi:1 0.1 371 /journal.pone.0097688.g009 



forms is more favorable than that of neutral P-CABs. Asp 139 in 
particular should be a very important binding site for protonated 
forms of P-CABs through hydrogen bonds and electrostatic 
interactions. The energy contributions of Asp 139 for protonated 
forms are all more favorable than those for neutral P-CABs. Thus, 
protonated form is the potent form of P-CABs. To enhance P- 
CABs pKa values and binding with the key residue Asp 139 would 
help increase the inhibition activity. The 100 ns disassociation 
molecular dynamics suggest that TAK-438 in protonated form 
could have the long dwell time and disassociate from H^, 
K'''-ATPase very slowly, which make it a likely competitor for 
PPIs. The findings in this work provide a better structural 
understanding of the binding sites in H"'",K^-ATPase and a basis 
for further rational design of novel P-CABs. 



Supporting Information 

Figure SI Ramachandran plot of the human H^,K^-ATPase 

model. 

(TIF) 

Figure S2 RMSD for the backbone atoms of the SCH28080 
and TAK-438-1 complexes in 100 ns disassociation molecular 
dynamics. 
(TIF) 

Figure S3 RMSF of each residue for SCH28080 and TAK-438- 
1 complexes in 100 ns disassociation molecular dynamics. 
(TIF) 

Table SI GHde docking Gscore, Glide energy and QM/MM 

energy (kcal/ mol) of P-CABs. 

(DOC) 




Figure 10. Interaction modes of SCH28080 (A) and TAK-438-1 (B) witli H%K^-ATPase after 100 ns disassociation molecular 
dynamics. Water molecules into the ion channel are shown in ball & stick model and the potassium ions near P-CABs are represented by blue ball. 
doi:1 0.1 371 /journal.pone.0097688.g01 0 
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Figure 11. The comparison of energy decomposition for residues 
disassociation molecular dynamics. 

doi:10.1371/journal.pone.0097688.g011 



in binding sites of SCH28080 and TAK-438-1 after 100 ns 
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